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We present numerical calculations of the impact of charge carriers-carriers interactions on the 
dynamical properties of quasiparticles such as renormalized velocity and quasiparticle inelastic scat- 
QQ ' tering lifetime in a gapped graphene sheet. Our formalism is based on the many-body GoW- 

, approximation for the self-energy. We present results for the many-body renormalized velocity 

■ suppression and the renormalization constant over a broad range of energy gap values. We find that 

' the renormalized velocity is almost independence of the carrier densities at large density regime. We 

also show that the quasiparticle inelastic scattering lifetime decreases by increasing the gap value. 
^ . Finally, we present results for the mean free path of charge carriers suppression over the energy gap 

' values. 

' PACS numbers: 71.10.Ay, 81.05.Uw, 71.45.Gm 
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I. INTRODUCTION 

The latest rival to succeed silicon's status is graphene, a single atomic layer of graphite make a truly tiny transistor 
^ ' to decrease the size and to improve the operational speed of the electronic devices. Silicon lost it's brilliant electronic 
^ properties in pieces smaller than about lOnm and practically the smallest silicon chips which has been used in silicon- 
based electronics is 45nm. Furthermore, silicon has some limitations in speed of operations. These restrictions lead 
to serious challenges for the Moore's law which states that the number of transistors can be placed inexpensively on 
an integrated circuit has increased exponentially, doubling approximately every two years. This growth cannot be 
(~| ' maintained forever and thus the search is on to find and use new materials which may be able to produce higher 
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performance and better functionality. 

The recent discovery^ of graphene in 2004, and its fabrication into a field-effect transistor-'^ , has opened up a new 
^-H • field of physics and offers exciting prospects for new electronic devices and apparently possible to come over those 
^ . aforementioned limitations. Graphene has instructive and unique physics with special intriguing electronic properties 
OO . which has attracted remarkable attentions.— First, the electronic properties of graphene are improved in sizes less 
, than lOnm . Second, the massless Dirac-like electrons move through graphene with almost near-ballistic transport 
(<~^ , behavior with less resistance because back-scattering is suppressed. Third, graphene is itself a good thermal conductor 
, such that graphene's thermal conductivity is about ~ 5.3 x lO'^ W/mK at room temperature which is greater than 
, the thermal conductivity of carbon nanotubes.'' fnterestingly, the mobility of carriers in graphene is quite high and 
J> ^ it is about fO^ cni^/Vs at room temperature."* It is important to note that the highest electron mobility recorded 
\ on the semiconductor junction H-Si(IIl)-vacuum FET is 8 x lO'^cm^/Vs at 4.2 K or the mobility of electrons in 
H ' junction Si-SiO2(100) MOSFET systems is 25 x lO'^cm'^/Vs at low temperature^, make graphene promising for different 
■ ■ ' applications in devices. 

Providing capability to control a type and density of charge carriers by gate voltage or by the chemical doping^ made 
graphene instructive for novel nano-electronic devices. However, a gapped semiconducting behavior would be more 
suitable for electronic applications. There have been some proposed in literature for a gap generation in graphene due 
to breaking of the sublattice symmetry by some substrates (such as SiC^, graphite^ and boron nitride^), to adsorbe 
some molecules (such as water, ammonia*^ and CrOa^*), spin-orbit interaction^^ and finite size effectJ^ In case, we 
are interested to carry out the microscopic theory to calculate some physical quantities of gapped graphene. 

Theoretical calculations of quasiparticle properties of electron in conventional two-dimensional electron liquid are 
performed within the framework of Landau's Fermi liquid theoryi^ whose key ingredient is the quasiparticle concept 
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and its interactions. As applied to the electron liquid model this entails the calculation of effective quasiparticle- 
quasiparticle interactions which enter the many-body formalism allowing the calculation of various physical prop- 
erties. A number of calculations considered different variants of the GoM^-approximation for the self-energy in 
two-dimensional electron ga a-'^^i-'^^i-'^'^^'-'^^i-'^^i^^i^"'^ from which density, spin-polarization, and temperature dependence of 
quasiparticle properties are obtained. 

There is a mechanism for quasiparticle scattering against quasiparticles because they interact through the Coulomb 
interaction. This is an inelastic process and induced a finite lifetime of the quasiparticles. The carrier lifetime in 
an epitaxial graphene layers grown on SiC wafers has been recently measured. Since experiments carried out their 
measurements on graphene placed on SiC, we expect that graphene was gapped. The experimental measurements are 
relevant for understanding carrier intraband and interband scattering mechanisms in graphene and their impact on 
electronic and optical devisesi^i^i 

In this paper we focus on the effect of energy gap on the renormalized velocity, the inelastic scattering lifetime 
of quasiparticles and the inelastic mean free path in gapped graphene sheets over the broad range of energy gap. 
Our formalism is based on the Landau- Fermi liquid theory incorporating the GqM^- approximation for the self-energy. 
These quantities are related to some important physical properties of both theoretical and practical applications such 
as the band structure of ARPES spectra^^, the energy dissipation rate of injected carriers^ and the width of the 
quasiparticle spectral functioni^ 

The contents of the paper are described briefly as follows. In Section II we discuss about our theoretical model 
which contains the effect of gap in the renormalized velocity of quasiparticles and the inelastic scattering lifetime Ti„ , 
of gapped graphene due to electron-electron interactions by using GoM^-approximation. Our numerical results are 
given in Section III. Finally, Section V contains the summery and conclusions. 

II. THEORETICAL MODEL 

Among the methods designed to deal with the intermediate correlation effects, of particular interest for its physical 
appeal and elegance is Landau's phenomenological theorj*i^ dealing with low-lying excitations in a Fermi-liquid. 
Landau called such single-particle excitations quasiparticles and postulated a one-to-one correspondence between 
them and the excited states of a non-interacting Fermi gas. He wrote the excitation energy of the Fermi-liquid in 
terms of the energies of the quasiparticles and of their effective interaction. The quasiparticle-quasiparticle interaction 
function can in turn be used to obtain various physical properties of the system and can be parameterized in terms 
of experimentally measurable data. In this paper, we will compute the energy gap dependence of the renormalized 
velocity, renormalization constant and the inelastic scattering lifetime of quasiparticle in a gapped graphene sheet. 

A. Quasiparticle renormalized velocity 

The dynamics of quasiparticles in a gapless graphene are described by two-dimensional (2D) massless Dirac Hamil- 
tonian H = hva ■ k, with eigenvalues Ssk ~ shvk, where s = +(— ) representing right- and left-handed helicity or 
chirality for the electrons and holes, respectively. Note that chirality is the same as helicity for the massless particles, 
t; = 10^ m/s is the Fermi velocity. As it has been shown before^^, contrary to conventional 2D electron systems, the 
interactions increase the velocity of quasiparticles in graphene because of interband exchange interactions and the 
difference between positive and negative energy branches due to the chirality. 

The dynamics of quasiparticles in a gapped graphene are described by 2D massive Dirac Hamiltonian given by 
H — hvpcr ■^+'mv^a3 with eigenvalues i?sk = sy^ {hvk)^ + where A = mw^ is the gap energy. Due to massive term 
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in the Hamiltonian, the chirahty differs from the hehcity and also the hehcity is conserved but is frame dependence. 
From the microscopic point of view, the quasiparticle energy can be calculated by solving the Dyson equation, 

fc?/ = 6k + 5Re[<5Er (k,^)]|_, (1) 

where 6k — ^-sk ~ ^'F is the energy of a quasiparticle relative to the Fermi energy. The Fermi wave vector in 
graphene is given by kp = {^T^n/ ggg^)^/"^ where gs — g^ — 2 are spin and valley degeneracy, respectively. The 
Fermi energy of gapless graphene is e-p — hvkp. The retarded self-energy of gapped graphene is Sj*^* and we define 
61^1'^'^ (k, Lu) — Sj'^*(k, w) — l]l'^'^{kp,0). In the on-shell approximation^^, on the other hand, the above equation must 
be solved by setting uj = 6k /ft- 

In the GqW- approximation^^, the self-energy of gapped graphene at finite temperature {(3 = l/{kBT)) is given by 



S^(k,^^«) - -^E / (^^^^'(k,k + q) (2) 



X 
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where the dynamic screened effective interaction is VF(q, i^m) — Vq/^il^ *^^m) and e{q, iflm) is the dynamical dielectric 
function and the bare Coulomb interaction is Vq = 2Tre'^ /nq where k is the averaged background dielectric constant of 
graphene is placed on a substrate. G^^\q,inm) — l/(*f^m — 6k/ft) is the standard noninteracting Green's function. 
The overlap function for gapped graphene F'"* (k, k -|- q) , is given by^S 

To evaluate of the zero temperature retarded self-energy, we decompose the self-energy into the line which is purely a 
real function and residue contributions, T,^'^* (k, uj) — T,'-^'^'^ (k, lu) + S^'^'* (k, w) , where S'*"^ is obtained by performing 
the analytic continuation before summing over the Matsubara frequencies, and S'''^'' is the correction which must be 
taken into account in the total self-energy, 



{k,u) ^ -^j j^VqF^^ (k,k + q) (4) 

dQ. 1 1 
2n e(q, iVl) oj + iO. — {k + (\) / 



and 



'(k,c.) = V / , (k,k + q 



oPq_ 

(2^)2e(q,c.-6'(k + q)A)^ 

X [e(c.-6'(k + q)A)-e(-6'(k + q)/ft)], (5) 

where the dynamic dielectric function is given by e(q, = 1 — VijX*-''^(g,ti;) in the random phase approximation 
(RPA) and x^'^\li^) is the noninteracting polarization function for gapped graphene. The noninteracting polarization 
function has been recently calculated on both along the imaginary and real frequency axis22i^. The noninteracting 
polarization function expressions along the real frequency axis^ are given in appendix A. 

Note that there are two independent parameters in the self-energy. One of them is the Fermi energy Ep, and the 
other is the dimensionless coupling constant agr = gggys'^ / kHv. The coupling constant in graphene depends only 
on the substrate dielectric constant while in the conventional 2D electron systems the coupling constant is density 
dependent. For graphene placed on SiC or graphite substrates, the coupling constant is about agr — 1. 
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The quasiparticle energy depends on the magnitude of k for the isotropic systems. Expanding Se'^^ to first order 
in fc — fcp, we obtain Se'^j^ ~ hv*{k — kp) which effectively defines the renormalized velocity as Hv* = dSe'^j^ /dk\k=kF ■ 
The renormalized velocity in the Dyson scheme is thus given by 

l-9^3?e[(5S';^*(k,c^)]|^=o,fc=fep ' ^' 

In the on-shell approximation, on the other hand, the renormalized velocity is given by v* /v = (1 + A^)~^/^ + 
v^^dk^e[SY^^_^*(]i,uj)]\^^Q^k=kF + (1 + A^)~^/^9aj3fJe[(5I]!j5*(k, w)]|^=o,fc=feF- The renormalized velocity in this approx- 
imation demonstrates qualitatively the same behavior obtained by the Dyson equation, Eq. ([5]) but its magnitude 
is larger than the one calculated within the Dyson schemci^ There is an ultraviolet divergence in the wave vector 
integrals of the line contribution in a continuum model formulated as discussed abovci^l We introduce an ultraviolet 
cutoff for the wave vector integrals, fee = Afcj? which is the order of the inverse lattice spacing and A is dimensionless 
quantity. For definiteness we take A = kc/kp to be such that 7r(Afci?)^ = {2n)'^/Ao, where = 3\/3ao/2 is the 
area of the unit cell in the honeycomb lattice, with uq ~ 1.42 A the carbon-carbon distance. With this choice, 
A ~ (.gn-i\/3/9.09)i/2 x 10^, where n is the electron density in units of 10 cm . 

An important quantity in the Fermi-liquid theory is the renormalization constant Z, defined as the square of the 
overlap between the state of the system after adding (or removing) of an electron with the Fermi wave vector and 
the ground-state of the system. The non-zero renormalization constant value is always smaller than the one for the 
normal Fermi-liquid systems and can be calculated explicitly as foUow^^ 

^ " 1 - d^ne[sj:-f{k,uj)]i=o,k=k, ■ 

We will show that Z is a finite number for gapped graphene and it confirms as well that the system is a Fermi-Liquid. 



B. Inelastic scattering lifetime 



In this subsection, we compute the inelastic scattering lifetime of quasiparticles due to carriers-carriers interactions 
at zero temperature and disorder-free for gapped graphene sheets. This is obtained through the imaginary part of 
the self-energy^ when the frequency evaluated at the on-shell energy. 

T-i(k) = r,„(k,e+k//i) = -^3mS-*(k,C+k/S), (8) 

where ri„(k, ^sk/^) is the quantum level broadening of the momentum with cigenstate |sk >. It is worthwhile to 
note that the expression of is identical with a result obtained by the Fermi's golden rule summing the scattering 
rate of electron and hole contributions at wave vector kj^ Note again that the total contribution of the imaginary 
part of the retarded self-energy comes from the residue term both intra- and interband contributions, 3r7iI]!J_^* (k, uj) = 
3mE[^(^^(k, Lj) + 3ml][^|g^(k, a;). However, the total contribution of the imaginary part of the retarded self-energy 
evaluated at the on-shell energy comes only from intraband term, 3mE!J_'^*(k, ^k/^) = S^'T-^rritraC'' Ck/^)- We will 
discuss about that with more details in the appendix B and C. 

We turn our attention to investigate the imaginary part of the retarded self-energy with more details. By starting 
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from Eq. ([5]), we end up to an expression for the imaginary part of self-energy which is given by, 



QmS7*(k,w) 



ires 
'intra 



(k, Uj) + ^mYi, 



^res 
'i7iter 



(k,c.) 



+ 



X 



X 




(9) 



where the imaginary part of the inverse dielectric function in RPA level is obtained by 



(10) 



[1 - Vg^ex^°^ (q, ujW + [Vg^-mx^"^ (q, w)]2 ' 



It is worth to note that the plasmon contributions in the imaginary part of self-energy comes from the zero-solutions 
of denominator in Eq. pO|) . 



We turn to a presentation of our main numerical results. We present some illustrative results for the quasipar- 
ticle dynamic properties such renormalized velocity, renormalization constant and inelastic scattering lifetime. All 
numerical data are calculated in the Dyson scheme at a = 1. 

The Fermi liquid phenomenology of Dirac electrons in gapless graphene^^i^l and conventional 2D electron liquid^S 
have the same structure, since both systems are isotropic and have a single circular Fermi surface. The strength of 
interaction effects in a conventional 2D electron liquid increases with decreasing carrier density. At low densities, the 
quasiparticle renormalization constant Z is small, the renormalized velocity is suppressed-^, the charge compressibility 
changes sign from positive to negative, and the spin-susceptibility is strongly enhanced'^^. These effects emerge from 
an interplay between exchange interactions and quantum fluctuations of charge and spin in the 2D electron liquid. 

In the 2D massless electron graphene, on the other hand, it has been shown^^i^li^^ that interaction effects also become 
noticeable with decreasing density, although more slowly, the quasiparticle renormalization constant, Z tends to larger 
values, that the renormalized velocity is enhanced rather than suppressed, and that the influence of interactions on the 
compressibility and the spin-susceptibility changes sign. These qualitative differences are due to exchange interactions 
between electrons near the Fermi surface and electrons in the negative energy sea and to intcrband contributions to 
Dirac electrons from charge and spin fluctuations. 

In this paper we have shown the results for gapped graphene which are determined values between the gapless 
graphene evaluated at A = and the conventional 2D electron liquid where A ^ oo. 

In Fig. 1, we have plotted the renormalized velocity as a function of carrier density for the various energy gap. As 
a result, we see that the impact of energy gap on quasiparticles velocity which is similar to the effect of impurity 
to that on graphene^. The renormalized velocity is almost density independent in gapped graphene at large carrier 
densities. The renormalized velocity reduces dramatically by increasing the energy gap especially in the low carrier 
densities. Importantly, the renormalized velocity becomes less than the bare velocity at large energy gap and low 
density values. It is physically accepted since the system tends to conventional 2D electron liquid by increasing the 
energy gap values. Note that in the conventional 2D electron systems, the renormalized velocity is suppressed by 
increasing the coupling constant or reducing the density. 



III. NUMERICAL RESULTS 
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We have shown the renormaHzation constant as a function of the energy gap in Fig. 2. The renormahzation 
constant enormously reduces by increasing the energy gap in mild densities, however it decreases quite slowly in high 
densities. 

Fig. 3(a) is shown the absolute value of 3mS'^*(k, uj) as from Eq. ([9]), evaluated at lu ~ ^k/^- By increasing the gap 
value, this function takes a finite jump at the wave number of the plasmon dip and at large A values, a discontinuity 
appears. The discontinuity is peculiar to 2D electron liquid.— It is absent in gapless graphene and starts to arise 
from the fact that the oscillator strength of the plasmon pole is non-zero at special k value for gapped graphene. 

Fig. 3(b) is clearly shown the behavior of the energy gap dependence of the inverse inelastic scattering lifetime. 
As it is argued in the Appendix B, the imaginary part of self-energy evaluated at the on-shell energy start from 
A — y/e^^+A^ and in case the results are truncated below that. The quasiparticle lifetime decreases by increasing 
the gap value and it is a clear difference between 2D massless Dirac electron and gapped graphene. Consequently, 
the inelastic scattering lifetime in graphene is always larger than the conventional 2D electron liquid. In the case of 
gapless graphene, scattering rate is a smooth function because of the absence of both plasmon emission and interband 
processesjSi nevertheless with generating a gap and increasing the amount of it, plasmon emission causes to arise a 
discontinuity in the scattering time, similar to conventional 2D electron liquidi^ We have thus two mechanisms for 
scattering of the quasiparticles. The excitation of electron-hole pairs which is dominant process at low wave vectors 
and the excitation of plasmon appears in a specific wave vector. We also see in Fig. 3(b) that the scattering rate is 
quite sensitive to the gap energy and the scattering rate increases by increasing the energy gap. 

In Fig. 4, we have depicted the inelastic mean free path Zi„(k) = t;*ri„(k), as a function of the on-shell energy 
for various gap energies. To this purpose we multiplied the results of Ti„(k) to a proper renormalized velocity. As a 
result the mean free path of a gapped graphene is shorter than that obtained for gapless graphene. Furthermore, the 
massless graphene has larger lin and it decreases by increasing the energy gap values. Note that the typical value of 
energy gap due to breaking sublattice symmetry is A = 10 — 100 meV corresponding the inelastic mean free path is 
lin = 20 — 50 nm which implies that the system remains in the semi-ballistic regimeiii2ii2iiiS. 

IV. SUMMERY AND CONCLUDING REMARKS 

In summary, we have studied the problem of the microscopic calculation of the quasiparticle self-energy and many- 
body renormalized velocity suppression over the energy gap in a gapped graphene. We have carried out calculations 
of both the real and the imaginary part of the quasiparticle self-energy within GoVF-approximation. We have also 
presented results for the renormalized velocity suppression and for the renormahzation constant over a wide range of 
energy gap. We have shown that the renormalized velocity for a gapped graphene is almost independent of the carrier 
density at high density. We have finally presented results for the quasiparticle inelastic scattering lifetime suppression 
over the energy gap and show that the mean free path of the charge carriers of a gapless graphene is larger than a 
gapped graphene one. In case the mean free path of charge carriers decrease by increasing the energy gap. 

A possible role of correlations including the charge-density fluctuations beyond the Random Phase Approximation, 
remains to be examined. 
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APPENDIX A: THE DYNAMIC POLARIZATION FUNCTION FOR A GAPPED GRAPHENE 



In this appendix we present the real and imaginary part of the noninteracting polarization function for a gapped 
graphene, which is calculated recently by Pyatkovskiy.'^" The dynamic polarization function for gapped graphene 
in the imaginary frequency axis is also calculated by us in Ref. [29|. Importantly, the noninteracting polarization 
function along the imaginary frequency axis can be obtained by performing analytical continuation from real axis and 
those results are the same<^ 

First, by introducing some following notations. 



hvk ' 



4A2 



(2 - a;o)cos"^(a;/a;o), 
(2 — Xq) cosh^"^(a;/a;o), 





— X? 


x"^ 


— Xq 


X^ 


Xq 



9± = 

Xo = 
G<(x) - 
G>(x) = 
Go(x) = 

the real part of noninteracting polarization function is given by. 



2lTVp 



and the imaginary part of noninteracting polarization function is given by, 

G>(g+)-G>(5-), 
G>(g+), 
0, 
0, 
0, 

-G<(5-), 
7r{2~xl), 



0, 




lA 


G<(5-), 




2A 


G<(5+) + 


G<g-), 


3A 


G<(.g-)- 


G<(5+), 


4A 


G>(.9+)- 


G>(5-), 


IB 


G>(5+), 




2B 


G>(5+)- 


G>(-5-), 


3B 


G>(-ff-) 


+ G>(5+), 


4B 


Go{g+) - 


Go 


5B 



f{k,uj) X < 



7r(2 
0, 



lA 
2A 
3A 
4A 
IB 
2B 
3B 
4B 
5B 



with the foUowings regions in the (fc, lu) space. 



lA 
2A 
3A 
4A 
IB 
2B 
SB 

4:B 

5B 



huj<EF~ y^h^V^jk - kpy + A2, 



y/h^vlik-kpy + A^l < huj < -Ef + ^Jh^v '^jk- 
hw < -Ef + v/^^w^(fc-fcF)2 + A2, 
-Ef + ^Jh?v'^{k + fcp)^ + A2 <huj < Hvk, 



A2, 



k < 2fcF, Vfi2w2fc2 + 4A2 <nuj <Ef + y/h^v^{k- kp)^ + A^, 
Ef + y/h^v^{k- kp^ + < hio < Ef + ^Jh^v^ jk + fcF)2 + A2, 

njjj>EF + y Jh'^v'^{k + kFY + A'^, 

k > 2kF, Vh^v^k^ + AA"^ < fk j<EF + y/h^ v^{k - fcp)^ + A2, 
hvk < huj < \/W-v'^k'^ + 4A2, 



(Al) 



(A2) 



(A3) 



(A4) 
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APPENDIX B: THE INTRABAND CONTRIBUTION OF SELF-ENERGY 

Since we are interested in quasiparticle properties, we therefore need only s = + contribution. Let us focus on the 
intraband contribution of the retarded self-energy. The second argument of the dielectric function in Eq. ([5]) ( by 
setting h = V — 1) is 



e+(k + q) =uj + Ef- v/fc2 + g2 + 2fc9Cos0 + A2 . (Bl) 



In this case, we change the variable (p and integrate it over y = y^k'^ + q'^ + 2kq cos + . Using the new variable, 
the intraband contribution of self-energy changes to 



27rKVfc2 + A2 Jo J^[k~qy+A^ ^Ak'^q^ - {y^ - - q^ - A^^ e{q,LU + EF-y) 

[e{Lu + Ef - y) - e{EF - y)]. (B2) 



We can now simplify the 0-functions further in Eq. (|B2p by considering the positive and negative regions of lo as 
follow 

1) uj +E-p — y > and Ep — y < : It implies that uj > , 

2) UJ +Ef — y < and E-p — y > : It implies that uj < Q . 

To consider the first case where a; > 0, the difference between the two 9-functions in Eq. (IB2p is equal to +1 if 



Ey <y<uj + EY and ^Jik-qf + A"^ <y< ^/{k + qf + ^ ■ (B3) 
Now we do need to find the overlap between these two intervals. We simply end up to inequivalent conditions which 



are q > k — y lu'^ + kp + 2u!y/kp + A"^, q < fc + y + fcp + 2ujy/kp + A^ and q > kp — k. Collecting everything 
together and using the fact that g > 0, wc finally find 



g2 rk+^uj'^+kl+luy/k^+A^ ^mm(w-(-^fc2-(-A2Y(fc+g)2+A2) 

,(k,w>0) = / dq / dy 

ZTTKVk^ + A^ Jmax{QMF-k,k-^ui^+k^+2uj^kl + A^) J max{^ k^ + A^ (k-q)^ + 

{VWr^ + yf-q' ^^^^ 



e(q, uj + Ep-y) ^ik^q^ ~ {y^ ~ - q^ - A^y 
By considering of the second case where w < 0, the difference between the two 0-functions in Eq. ^ is equal to —1 if 



Ef +UJ <y < Ef and (k - qf + A^ <y < y^^k + qf + A? (B5) 
As what we did before, we calculate overlap between intervals and thus we find q > k — kp and q < k + k-p, 



q > y ui'^ + kp + 2Ljy^kp + A^ — k. Putting everything together and using the fact that q > we finally find 
^:tra{K A~Ep<u;<0) = — 7^5== /'^'' ^ —= dq 

27rKVfc + A^ Jmax{a,k-kF,^uj^+k^+2u:y/k^+A2~k) 



|.mi«(^fe|,+A2,^(fe+g)2+A2) 

X / dy 

Jmax{0,LJ+y'kl+A^,y'(k~q)^+A^) 

e(q, uj + Ef- y)y^ik^q^ - {y^ - k^ - q^ - A2)2 
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r-fc+fcF /.mm(^fe|+A2,^(fe+g)2+A2) 

S[-,,(k,a; < -A - Ep) = == / dq / dy 



(B7) 



e(q, lu + Ef~ y) v/4fc2g2 - (y2 _ /,2 _ ^2 _ ^2)2 
The real and imaginary part of intraband contributions can be computed. 

APPENDIX C: THE INTERBAND CONTRIBUTION OF SELF-ENERGY 

Now, we focus on the interband contribution of the retarded self-energy. The second argument of the dielectric 
function in Eq. ([5]) is 



w-^_(k + q) ^ uj + Ef + + + 2kqcos(l) + . (CI) 



We change variable y = y^k^ + q'^ + 2kq cos cj) + A^, then we find 



yres . f^^ , f^/i^^^ dy q' - {y - ^/FTA^) 



27rKVl2~:rA2 7o 7^ (fc^g)2+A^ ^4fc2g2 „ (y2 _ fc2 „ ^2 _ ^2)2 e(q, + + y) 

X [e(cj + £;F + y) -1]. (C2) 

Note that '^inter '^^^ t)e non-zero if w + E'p + t/ < and y > 0- It means that lu < —Ep. In this case the difference 
between the two 0-functions in Eq. ([5]) becomes -1 if 



0<y<-{u; + EF) and y^{k - qf + A^ < y < y^ik + qf + A"^ . (C3) 
Now we do need to find the overlap between these two intervals. We end up to inequivalent conditions that q > 



fc — y Lj2 _|_ fc2 _|_ 2ujy^kp + A2 and q < k + y uj'^ + kp + 2uj kp + A"^ . Putting everything together and using the fact 
that g > we finally get 



e 



2 



fc+-y/w2 + fe|+2w^fc|+A^ <.min(^(fe+g)2+A2,-(w+^fe| + A2)) 



E-L,(k,c. < -Ef) = — yf^=^ / , —= dq / dy 

A^v r i-i ^ max (0,k-^iJ^+kl,+2iUy/kl + A-2) J^(fc-q)2+A2 

g2_(y_^/^J-^)2 ^^^^ 



e(q, -I- ep + y)^j4.k'^q'^ - (y2 ^ fc2 - g2 _ 



If we want to calculate 3TOS]!j_'^*(k, ^_|_(k)) needed for computing the quasiparticle lifetime, wc will only need the 
intraband contribution of the self-energy since the interband contribution is zero. 
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FIG. 1: ( Color online) (a): The renormalized velocity as a function of density for the various energy gaps at agr = 1. (b); 
The renormalized velocity as a function of the energy gap for the various densities. 
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FIG. 2: ( Color online) The renormalization constant as a function of the energy gap for the various densities. 
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FIG. 3: ( Color online) (a): The absolute value of the imaginary part of the retarded self-energy on the energy shell as a 
function of the wavevector for the various energy gaps; (b): The inelastic quasiparticle lifetime(ri„) in graphene as a function 
of the on-shell energy for the various energy gaps at n = 5 X lO'^^cm"^. 




FIG. 4: ( Color online) The quasiparticle mean free path as a function of the on-shell energy for the various energy gaps at 
n = 5 X 10"cm-^ 



